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Abstract: We present the temporal changes of the diffusion coefficient K of galactic cosmic rays (OCRs) at 
the Earth orbit calculated based on the experimental data using two different methods. The hrst approach is 
based on the Parker convection-diffusion approximation of GCR modulation IIl:i.e. K ^ Vr/dl where dl is the 
variation of the GCR intensity measured by neutron monitors (NM), V is the solar wind velocity and r is the radial 
distance. The second approach is based on the interplanetary magnetic field (IMF) data. It was suggested that 

5 

parallel mean free path Aj| can be expressed in terms of B as follows 0-11: A||- where SB is the standard 
deviation. In our calculations we used an approximately equivalent expression A|| ■^. Using data of the product 

A|| V^n of the parallel mean free path A[| and radial gradient V^n of GCR calculated based on the GCR anisotropy 
data (Ahluwalia et ah, this conference ICRC 2013, poster ID; 487 Q), we estimate the temporal changes of 
the radial gradient V^n of GCR at the Earth orbit. We show that the radial gradient V^n exhibits a strong solar 
cycle dependence (11-year variation) and a weak solar magnetic cycle dependence (22-year variation), being in 
agreement with the previous other calculations and with PIONEER/VOYAGER observations. 
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1 Introduction 

The galactic cosmic rays (GCRs) transport in the helio¬ 
sphere is governed by the four important processes: out¬ 
ward convection by the solar wind, inward diffusion, parti¬ 
cle drifts (gradient, curvature and on the neutral sheet) in 
the turbulent interplanetary magnetic field (IMF) and adia¬ 
batic cooling. Estimation of the local electromagnetic con¬ 
ditions near the Earth orbit is possible by establishing mod¬ 
ulation parameters as diffusion coefficient, density gradi¬ 
ents etc. It is especially essential when in situ measure¬ 
ments are absent. It is the basic knowledge needed to study 
space weather prediction. A major issue in GCR transport 
research and space weather studies is how GCR particles 
propagate through the heliosphere, and how interact with 
the interplanetary space especially in the inner heliosphere 
near the Earth’s orbit. The essential significance of charac¬ 
terizing GCR propagation is evident, because this will lead 
to a practical capability in space weather forecasting which 
has important consequences for life and technology on the 
Earth and also in the interplanetary space. 

The hrst theoretical description of cosmic ray transport 
coefficients was done by Jokipii ||6l by formulation the 
quasilinear theory (QLT) for GCR diffusion. One limiting 
assumption of the QLT is that in the guiding center approx¬ 
imation transport of the GCR particles is not perturbed by 
the IMF turbulence. This assumption is inaccurate especial¬ 
ly for the highly anisotropic strong turbulent heliosphere. 
This classical approach has been improved by higher-order 
theories of GCR particles turbulent how. Thus the theo¬ 
ries considering nonlinear effects have been introduced Q- 
IfTTl . A validity of the QLT for the GCR particles of the en¬ 
ergy > \GeY is conhrmed by the weakly nonlinear theory 
(WNLT) IfTOl . nonlinear parallel diffusion theory (NLPA) 
lfT2l and in papers lfT3l fT4l fTTl (see e.g. GS). 

However, for selecting correct set of modulation param¬ 
eters used in theoretical modelling (especially diffusion co¬ 
efficients, etc.), one criterion remains the most important. 


if it is possible, to estimate them from the experimental da¬ 
ta, supported by appropriate theory. 

Observations of GCR intensity and anisotropy by neu¬ 
tron monitors (NMs) and IMF Huctuations can be success¬ 
fully used for establishment of various parameters charac¬ 
terizing modulation of GCR by the solar wind. In this pa¬ 
per we present the temporal evaluation of the parallel dif¬ 
fusion coefficient of GCR particles for rigidities to which 
NMs respond. Parallel diffusion coefficient equivalent 

to parallel mean free path (MFP) A|| = determined by 
physical properties of interplanetary medium, is a very im¬ 
portant parameter to study the transport of energetic parti¬ 
cles in the heliosphere, especially for a solar event (SEP) 
connected with the space weather prediction HSlIIll. 

Using data of the product A|| V^n of the parallel MFP A[| 
and radial gradient V^n of GCR calculated based on the 
GCR anisotropy data (Ahluwalia et al., this conference I- 
CRC 2013, poster ID: 487 0), we estimate the temporal 
changes of the radial gradient V^n of GCR at the Earth’s or¬ 
bit. As a hnal point, determination of the parallel diffusion 
coefficient (equivalent to parallel MFP A||) of GCR par¬ 
ticles according to Quenby ifTSi and Hedgecock ifT^ for¬ 
mulas will be performed for minimum conditions of solar 
activity. 

2 Convection-diffusion approximation 

It has been shown that ^ 75 — 80% of the 11-year varia¬ 
tion of the GCR intensity can be interpreted based on the 
diffusion—convection model of GCR propagation ll20ll2TII . 
So, the long-term variations of GCR intensity can be de¬ 
scribed by the Parker transport equation, invoking the 
isotropic convection-diffusion approximation fT2\ . In s- 
cope of this approximation one can calculate changes of 
the parallel diffusion coefficient Kn as follows: 
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Fig. 1; Temporal changes of the annual SSN, GCR intensi¬ 
ty for Moscow NM, solar wind velocity V, magnitude B of 
the IMF, product VB and estimated parallel diffusion coef- 
hcient/fji for 1965-2011. 
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Fig. 2; Temporal changes of the 13-month smoothed mag¬ 
nitude B and Bx, By and components of the IMF for 
1965-2011. 
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Fig. 3: Temporal changes of the 13-month smoothed SSN 
and parallel MFP A|| calculated according to expressions 
(3) and (4) for 1965-2011. 


where dl is variation of the GCR intensity, V solar wind 
velocity and r radial distance. 

Figure [T] shows a plot of the temporal changes of the an¬ 
nual sunspot numbers (SSN), GCR intensity for Moscow 
NM, the solar wind velocity V, the magnitude B of IMF, 
product VB and the estimated parallel diffusion coefficient 
/Til according to expression (1) for the period of 1965-2011. 
One notes that the parallel diffusion coefficient exhibit- 
s 11-year variation, but a stronger solar polarity depen¬ 
dence (~ 22-year variation); a significant increase is ob¬ 
served in the minimum epochs of solar activity, especially 
in the A < 0 magnetic polarity period. An anomalous in¬ 
crease of 7^11 for the recent solar minimum 23/24 is clearly 
seen, as well. 

3 Parallel mean free path 

It has been suggested that parallel mean free path (MFP) 
A[| can be expressed in terms of interplanetary magnetic 
held B as follows El-ll: 



( 2 ) 


where dB is the standard deviation. In our calculations the 
formula (2) is replaced by an equivalent expression; 


A||OC 


B 
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We calculate parallel MFP A|| from the 27—day running av¬ 
erages of the observed radial Bx, azimuthal By, and latitu¬ 
dinal Bj components of IMF according to the formulas: 


Bx ^ By B^ 
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(4) 


Data sets of the IMF magnitude B and Bx, By, B^ com- 
ponets used in this calculation are presented in hgure |2l 
The results of MFP calculations for 13-month smoothed 
data for the corresponding components and magnitude of 
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Fig. 4; Time variation of the annual GCR radial gradient Fig. 5; The time variation of the resonant frequency / for 
V^n for 1965-2011 with errors bars. the 1965-2011 period. 


the IMF according to expressions (3) and (4) are present¬ 
ed in figure |3] One notes that the MFP oscillates with a pe¬ 
riod of ~ 11 year solar activity cycle with a significant in¬ 
crease in the minimum periods of solar activity. Also, MFP 
is strongly polarity dependent in accord with the drift theo¬ 
ry with a considerable enhancement especially in the mini¬ 
mum epoch of solar activity in the A < 0 magnetic polarity 
period. 
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4 Radial gradient of GCR 

On the basis of the long term changes of the product A|| V^n 
calculated based on the GCR anisotropy data (Ahluwalia 
et al., this conference ICRC 2013, poster ID; 487 Q) and 
parallel diffusion coefficient of GCR found above (fig¬ 
ure [T]), we estimate also the radial gradient V^n of GCR 
at the Earth’s orbit. The results of our calculations are pre¬ 
sented in figurelH One notes that the radial gradient V^n of 
GCR oscillates with a period ^ 1 lyear solar activity cycle 
with a weaker solar polarity dependence, being in agree¬ 
ment with the previous calculation reported by Chen and 
Bieber ||23]| and with PIONEERA'OYAGER observations 
El. 


5 Transport coefficient for minimum 
conditions of solar activity 

Transport coefficients (e. g. diffusion coefficient of GCR) 
may be derived from a precise knowledge of the regular in¬ 
terplanetary magnetic field values and its fluctuations (tur¬ 
bulence) M . The derivation of the parallel diffusion coef¬ 
ficient given by Jokipii lH |25]| and Hasselmann and Wib- 
berenz ESj is best illustrated in a simple way by following 
the Kennel and Petschek EtII formulation given in the con¬ 
text of magnetospheric particle scattering. 

In order to accurately obtain the parallel mean free path 
(diffusion coefficient) of GCR particles in the heliosphere, 
a method of power spectrum density of the interplanetary 
magnetic held turbulence has been used im. In this paper 
we compare the values of parallel MEP obtained by means 
of formulation of Hedgecock 1(191 and Quenby ifTSl . We 
consider frequency range 10^® — 10^^Hz, responding for 
modulation of the GCR particles detected by NMs. 

The appropriate parallel MEP can be expressed in terms 
of the power spectrum density of the interplanetary mag¬ 
netic held fluctuations according to Hedgecock 1(191 and 
Quenby respectively; 


Fig. 6; Values of the parallel mean free path A|[ calculated 
based on expression (5) according to Hedgecock HU for 
the consecutive minimum epochs of solar activity with 
different signs of global magnetic polarity for the period 
of 1975-2011. Each point corresponds to the average value 
for three years around each solar minimum period. 


, 2v{v + 2)cR^ 

II 9V-P{f) 


(5) 


, VB^ 1 

Where V is the solar wind velocity, B magnitude of IME, R 
- magnetic rigidity of GCR particles to which NM respond 
(in this case R = 15GV), / is the resonant frequency of 
GCR scattering, P{f) is the power spectrum density at 
the resonant frequency / with the spectral index v for the 
frequency range 10^® — The time variation of the 

resonant frequency / = for the period 1975-2011 is 

shown in figure |5] Figure (5] shows clear 11 —year variation 
in accordance with the solar activity cycle. 

Figures 12 and |7] present the values of the parallel mean 
free path A|| calculated based on expression (5) and (6) ac¬ 
cording to Hedgecock 1(191 and Quenby l(T8l . respectively, 
for the consecutive minimum epochs of solar activity with 
different signs of global magnetic polarity for the period 
of 1975-2011. The parallel mean free path A|| is calculat¬ 
ed based on the transverse components By and B^ of the 
IMF. Each point corresponds to the average value for three 
years around each solar minimum period (e.g., 1986 cor¬ 
responds to the time interval 1985-1987). Calculations for 
parallel MFP according to both formulas (5) and (6) are in 
good agreement with each other. One can see that parallel 
MFP calculated based on the By transverse component of 
the IMF is strong polarity dependent with large increase in 
A < 0 for 1985-1987. Unfortunately this statement for By 
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Fig. 7; Values of the parallel mean free path A|| calculat¬ 
ed based on expression (6) according to Quenby im for 
the consecutive minimum epochs of solar activity with d- 
ifferent signs of global magnetic polarity for 1975-2011. 
Each point corresponds to the average value for three years 
around each solar minimum period. 

component is not satished in the last minimum 2007-2009 
(A < 0). On the other hand in the last minimum period with 
record level of the OCR intensity ever measured by NMs, 
calculations for transverse component show an increase 
in the changes of the parallel MFR 

6 Conclusions 

1. The parallel diffusion coefficient , obtained based 
on the isotropic convection-diffusion OCR modula¬ 
tion model, generally displays ^ 11-year variation, 
but with strong polarity dependence 22 years). A 
signihcant increase is observed in the minimum e- 
pochs of solar activity, especially in the A < 0 mag¬ 
netic polarity period. An anomaly increase of /Tn in 
recent solar minimum 23/24 is clearly seen, as well. 

2. We calculate parallel mean free path A|| of GCR 
based on the experimental data of the IMF. Its value 
is polarity dependent in accord with drift theory and 
oscillate with a period of ~ 11-years solar activity 
cycle. 

3. On the basis of the long term changes of the GCR 
anisotropy we show the ^ 11-year variation of the 
radial gradient V^n of GCR being in good agreement 
with the PIONEER/VOYAGER observations. 

4. Parallel mean free path A[| calculated based on the 
By transverse component of the IMF is strong po¬ 
larity dependent with large increase in A < 0 for 
1985-1987. In the last minimum epoch of solar activ¬ 
ity with record level of the GCR intensity ever mea¬ 
sured by NMs, calculations for B^ transverse compo¬ 
nent of the IMF show an increase in the changes of 
the parallel mean free path. 
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